PAMELA through a Magnetic Lense 

J. p. Robert^] 
CCPP, 4 Washington Place, NYU, New York 



O Abstract 

(N 

^^ The PAMELA satellite has observed an excess of positrons over electrons in the 

^ energy range 1 — 100 GeV that increases with energy. We propose that the excess 

" is not due to a change in the local interstellar spectrum, but is due to heliospheric 

'^ modulation. We motivate this from the known form of the heliospheric magnetic 

field and predict that the excess will disappear when we enter a period of solar 
I""! maximum activity. 

m 
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Ph 1 Introduction 

6 

in The PAMELA satellite [I] has observed a rise in the positron to electron fraction that 

i__i increases with energy above 10 GeV . This has been interpreted as proof of a rise in the 

positron fraction in the local interstellar spectrum of cosmic rays. Such a result requires 

^ that there be some hardening of the positron spectrum at these energies and has caused 

^ a great deal of discussion of possible sources from decaying or annihilating dark matter 

\^ to acceleration from a population of pulsars. 

^. 

\f^ We propose instead that this rising fraction is due to the heliospheric magnetic field. 

O The heliospheric magnetic field is ordered over large distances during solar minimum p] [3j 

and is well described up to medium latitudes by a Parker spiral field [Ij. During solar 
minimum this means that for a particle to propagate directly to earth it must traverse a 
large scale ordered magnetic field that is primarily oriented perpendicularly to the radial. 
Direct transport is impossible for low energy particles [E < \ GeV^ as their Larmor 
^ radius is small even in the outer heliosphere. They instead reach the central solar system 

primarily by following the magnetic field lines along with a complex process involving 
convection, diffusion, drifts and energy loss. At high energies [E > \ TeV), particles 
can traverse the ordered magnetic field with little deviation. At these energies we see the 
interstellar spectrum directly. 

If the magnetic field were the only effect on the particle transport, and if the magnetic 
field were time invariant, then the switch from one set of particles paths to the other 
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would have no effect on the ffux by Liouville's theorem. However the magnetic field is 
time variant. The magnetic field lines are carried by the supersonic solar wind and travel 
out towards the edge of the heliosphere at 400 km/s. This causes large scale convection of 
cosmic rays that results in a deficit in the cosmic ray density in the central solar system 
when compared to the interstellar spectrum. This has been well measured for cosmic rays 
with rigidities up to 2.5GV |3],|5]. At these rigidities, the particles generally have a very 
long path length to the center of the solar system as they follow the magnetic field lines. 
Therefore they have a low radial velocity, and are greatly effected by the convection of 
the solar wind. The solar magnetic field also varies on longer timescales, affected by the 
27 day rotation period of the sun, and by the reversal of the solar magnetic field every 11 
years. 

At higher rigidities particles can travel directly to the central solar system at the speed 
of light. Their path length within the heliosphere is short and they are almost unaffected 
by convection from the solar wind. Therefore at high energies we expect there to be 
no radial gradient of cosmic rays as they should free stream through the solar system. 
From these two observations - the two separate transport paths for low and high energy 
particles, and radial gradient of low rigidity cosmic rays from the convection of the solar 
wind - we can already say that there must be a hardening of the cosmic ray spectrum 
between rigidities of 1 GV and 1000 GV. This will be due to the change over from a 
locally suppressed cosmic ray flux to the interstellar flux. 

To explain the PAMELA signal we need one further element. We require that the 
change over from the local flux to the interstellar flux be charge asymmetric. We already 
know that the orientation of the magnetic field results in asymmetric cosmic ray transport. 
Indeed this is the cause of the low energy heliospheric modulation below a few GeV^[L6j . We 
argue that this is also true for high energy particles and that the asymmetry should become 
more pronounced with increasing energy. We use a simple model of the heliospheric 



magnetic field to show this in section 4.1 This exhibits the primary features of the 



asymmetry that we are interested in. We then extend the study to consider the asymmetry 



in a more complex realisation of the field in section 4.2 



By postulating a heliospheric cause of the positron excess we can immediately make 
a clear prediction. We predict that the excess should disappear when the solar system 
enters a period of solar maximum activity which it is due to do shortly. In addition, in 
the opposite solar cycle when the sun's magnetic field is oppositely aligned, we would 
expect to see the opposite modulation - giving a rising excess of electrons over positrons 
in the energy range 1 — 100 GeV. We will also show that the tilt angle of the heliospheric 
current sheet (HCS) is an important feature of the model, and predict that any temporal 
variation of the latitudinal extent of the current sheet should be matched by a variation 
in the observed positron excess, with a time delay due to the time taken for magnetic field 
features at the sun to propagate out to the heliopause in the solar wind. These predictions 
are robust features of our model and testable in the near future. 



2 The heliospheric magnetic field 

In 1958 E. N. Parker proposed that the solar system magnetic field took the form of an 
Archimedes spiral [3]. The sun is a rotating dipole and the magnetic field is carried out 
into the solar system by the supersonic solar wind. As the solar wind travels out from 
the sun its rotational velocity slows with respect to the plasma at the surface of the sun, 
causing it to be retarded and dragging the magnetic field into a spiral. This continues 
out to the termination shock at around 80 AU where the solar wind becomes subsonic. 



The heliospheric magnetic field is described in the Parker model by: 
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where {r,6,(j)) are heliocentric spherical coordinates, Bq is the magnetic field at tq, V is 
the solar wind speed, b is the radius at which the field is purely radial and ^^0 is the 
equatorial rotation of the sun. 

The Parker model is very successful at modeling the magnetic field at moderate lat- 
itudes during times of solar minimum activity. More sophisticated models that perform 
better at high latitudes and times of greater solar activity have since been proposed. One 
widely used model is due to Fisk [7]: 
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where ug and u^f, are the differential rotation rates in the 6 and (f) direction respectively 



In both cases the magnetic field is essentially a spiral. The field lines are oriented 
radially close to the sun's surface. At large radii they are predominantly azimuthal. The 
winding of the magnetic field causes it to remain strong to large radii, with the strength 
of the azimuthal component falling off as 1/r. The solar magnetic field is also oppositely 
aligned in the northern and southern hemisphere, with a sheet of zero magnetic field 
separating the two. This is known as the heliospheric current sheet (HCS). 

The sun's magnetic poles are not aligned with its axis of rotation. This means that 
the HCS is tilted at the sun and that this tilt varies as the sun rotates. These variations 
are carried out into the solar system by the solar wind they cause the current sheet to 
be rippled. As the variations are carried out radially the latitudinal extent of the current 
sheet is described by a tilt angle dt- At times of solar minimum, the tilt angle is small, 
often around 5 — lOo. At times of solar maximum the tilt angle grows to large values, 
before the solar polarity flips and the new polarity is carried out to the heliopause by the 
solar wind [9]. 



The Magnetic Field in the Hehosphere 
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Figure 1: A schematic of the form of the heliospheric m,agnetic field in the current solar 
minimum. 

We are currently in a period of solar minimum and the magnetic field is shown schemat- 
ically in Fig. [l| The heliospheric magnetic field has been well mapped by a number of 
satellites as well as by ground based observatories. The current field configuration means 
that the field lines are directed in towards the sun in the northern hemisphere and outward 
from the sun in the southern hemisphere. The local magnetic field at earth in around 
3.3 nT. For a review see [3]. 



3 Cosmic ray propagation through the hehosphere 



To consider the motion of galactic cosmic rays through the hehosphere let's first state the 

equation for the Larmor radius of a relativistic particle in a magnetic field in convenient 

units: 

_ 2.2 X 10-V(GeV/c) 

'^^^^^^" \q\B{nT) ' ^^^ 

If we take the magnetic field strength (at low latitudes) to be \B\ = Be/t{AU) and 
the local field strength at earth over the course of the PAMELA data to be Be = 3.3 nT 
then we can approximate the Larmor radius of a particle at low latitudes to be: 



rg{AU) = 6.67 x 10-^p^{GeV/c)r{AU), 
where r is just the distance from the center of the solar system. 



(4) 



From this we can see that a singly charged particle with a rigidity of 1 GV has a 
Larmor radius of less that 1 AU even in the weak magnetic field near the termination 
shock at 80 AU. There is no means for such a particle to propagate directly to the center 
of the solar system. At the other extreme, a 1 TeV particle has Vg = 6.7r at all r. As 
Tg ^ r at all points on its path, it can penetrate the center of the solar system with 
relatively little deflection in its trajectory (around 10°). At energies above a few tens 
of TeV the deflection becomes negligible and we can truly say that any effect of the 
heliospheric magnetic field should disappear. 

In between these two limits we must consider two separate populations of cosmic rays: 

1. those that travel to the center of the solar system by traveling along magnetic field 
lines 

2. those that can penetrate the center of the solar system by traveling perpendicular 
to the magnetic field. 

At high energies we know that the flux of cosmic ray must be approximately homo- 
geneous and isotropic. The flux that we see should correspond directly to the interstellar 
spectrum. 

The same cannot be said of the low energy cosmic rays. The transport of suh-GeV 
cosmic rays within the solar system is described by the transport equation - flrst proposed 
by Parker in 1965 |10j : 

f)TT 1 f) 

_ = -V ■ (CVU) - V ■ ((v,)f/) + V ■ (/.(^^ ■ \/U) - 3 ^(PV ■ Vf/), (5) 

where U is the density of cosmic rays as a function of position, momentum and time,V 

represents the solar wind velocity, C = 1 — {l/3U)d/dp{pU) is the Compton-Getting 

coefficient, (v^) is the averaged drift velocity and k^^^ is the symmetric part of the diffusion 

tensor: 

'k|| 

K = \ K± Kt \ ■ (6) 
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The terms on the right hand side of the transport equation represent the change in 
particle density over time. The flrst term represents convection caused by the solar wind 
carrying particles out of the solar system towards the edges. This effect is large when 
the average radial velocity of a particle is of the same order or smaller than the solar 
wind speed. The second term quantifles the effects of drifts - where a variation in the 
magnetic fleld strength allows for a motion of the particles in the direction of V x B. These 
drifts are large scale particle motions and are modeled as effective sources or sinks in the 
system. They only affect the particle density if VU is non-zero. Convection ensures a 



radial gradient in the cosmic ray density and makes drifts a critical element of low energy 
cosmic ray transport. The third term covers diffusion - modeling the scattering of particles 
off inhomogeneities in the magnetic field. The terms in the diffusion matrix are not well 
known. For low energy transport, perpendicular diffusion is taken to be significantly less 
efficient than parallel diffusion [k± ~ 0.02fi;||). The final term represents adiabatic energy 
changes. The solution of this equation in 3D requires sophisticated numerical simulations 
which we will not discuss further here. For a review of recent work in this area see p,[TT]. 

For our purposes it is enough to note that the transport of low energy cosmic rays is 
very different from the transport of high energy cosmic rays. The average radial velocity 
of low energy cosmic rays is small - through perpendicular diffusion and drifts - and 
therefore the convection effects of the solar wind create a radial gradient in the cosmic 
ray density. This is well measured for low energy cosmic rays {E < 200 Me\^)and was 
analysed over different solar minima and maxima in |12]. The gradient has also been 
measured for protons and electrons with higher rigidities along the Ulysses trajectory [5]. 
In the inner heliosphere electrons and protons with rigidities of 2.5GV were shown to have 
a radial gradient of 2 — 3%/AU at solar minimum. In contrast we know that there cannot 
be any radial gradient at very high energies when the cosmic rays free stream through the 
solar system. 

We will take the approximation that the local density of cosmic rays that travel to the 
center of the solar system along the field lines to be suppressed with respect to the total 
interstellar flux by a factor A: 

FiAu = AFsoAu (7) 

This remains true up to 450 GeV when the particles can penetrate the central solar system 
directly. Clearly this is an approximation as particles will penetrate further into the 
solar system as their energy increases, affecting the overall flux before we reach 450 GeV. 
However within the energy we are interested in for PAMELA, we expect the approximation 
to be good as the Larmor radii remain small compared to the distance travelled. 

This energy dependent flux deficit would have little bearing on the PAMELA result 
if there was a smooth transition from one propagation method to the other. However, 
the form of the magnetic field allows some particles to penetrate directly into the central 
solar system even at low energies, and this effect is charge asymmetric. 



4 The magnetic lense 

4.1 A flat current sheet 

Let us consider a magnetic field that is purely azimuthal and has a strength that falls off 
as Be/t{AU). Let us further take the approximation of a fiat current sheet in the ecliptic 
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Particle Paths in a Simplified Heliosplieric Magnetic Field 




Figure 2: We show paths of positively charged particles that cross the ecliptic plane in the 
approximation of a fiat current sheet and a simplified magnetic field. Note that it does not 
matter at what angle the particle crosses the ecliptic plane, the overall motion is always towards 
the center of the solar system. 



plane. In this magnetic field, the Larmor radius of the particle is given by: 

2.2 X 10-V(G'ey/c)r(v4f/) 
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In Fig. [2] we show a schematic of the paths of positively charged particles that cross 
the ecliptic in the current magnetic field configuration. For any particle that crosses the 
ecliptic, its overall motion is towards the center of the solar system. For particles of the 
opposite sign the motion of a particle that crosses the ecliptic is away from the center of 
the solar system^ 



^We are certainly not the first to notice this effect. Indeed it was first shown to exist by Levy in 1976 
[131 [13] . He noted that there should be a large charge asymmetry and that this should exist up to large 



As the Larmor radius decreases with the distance from the center of the solar system 
we can define the volume within which the particles can travel directly to the center of 
the solar system. The fractional volume is just given by the opening angle: 

r \q\BE{nT) 

where / = sin 6 gives the fraction of the total solid angle that allows particles to travel 
through the ecliptic in this manneiQ 

Clearly this formula is not valid at large angles as we would need to carefully treat the 
varying Larmor radius of the particle over the course of one rotation due to the latitudinal 
variation in the magnetic field. Here we are considering only particles within the ranges of 
the PAMELA measurements, which corresponds to an opening angle of 33 degrees, small 
enough that these corrections can be neglected in our rough approximation. 

Within this volume there is a population of cosmic rays that arrive at the center of the 
solar system traveling perpendicular to the magnetic field lines rather than along them. 
This population has a radial velocity close to the speed of light and thus the effects of 
convection should be negligible. In addition, the overall flow of the opposite sign particles 
out of the solar system results in a relative decrease in the number density of oppositely 
charged particles. We model the result of this by taking the fluxes of the two populations 
of cosmic rays to be: 

F,Au{e^) = {A + f)F80Au{e^) 

F,Au{e-) = AFsoAuie-) (10) 

where / is the fraction of the particles in the heliosphere on a trajectory that allows them 
to travel across the HCS, and A is the suppression of the local flux with respect to the 
interstellar flux. El 



rigidities - in his case he took it to be up to 30GV, with the high energy cut-off being due to reduced 
radial gradients at these energies. Fohowing this work this effect was considered using the formahsm of 
drifts. At low energies the volume of particles that can drift in or out of the central solar system along 
the current sheet is small. Instead the drift of particles down from the poles dominates, resulting in the 
excess of electrons over positrons at low energies. It was argued that the oscillation of the current sheet 
reduced the effectiveness of the drift along the current sheet, resulting in no modulation above a few 
GeV. We will address this point in section 4.2 

^As we saw in Fig. I2I particles that are up to 2rg from the ecliptic can cross the ecliptic and drift 
towards the center of the solar system. However as you consider particles at increasing distances from the 
ecliptic, the range of angles that allow the particle to reach the ecliptic decreases. The fraction of particles 
at a height d from the ecliptic that have trajectories that can cross it is given by l/7rcos~^(d — Tg). When 
we integrate over this distribution we find that 1/2 of all particles in the volume defined by the angle 
—26 — >■ 29 have trajectories that can cross the plane. This gives the total fraction of the particles to be 
sin 9. 

■^Here we consider electrons and positrons, but we note that the same arguments can be made for 
other cosmic ray species. Specifically we would expect to see a deficit of anti-protons in the current solar 
minimum. For protons and anti-protons we cannot approximate rigidity with energy and we expect the 



Fit to PAMELA data with a flat current sheet 
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Figure 3: Using Eq. \1C\ we fit the PAMELA spectrum. We use a standard background spectra 
generated using Galprop \15\ Ii6]/ . 



Using this model we fit to the PAMELA data in Fig. [3j This fit requires that the local 
flux be 23% of the interstellar flux, though the large error bars on the high energy points 
allow for a wide range of values (14% to 30%) for the local suppression. In addition there 
is substantial uncertainty in the background spectrum which will have an effect on any 
flt. Nevertheless, we note that though the flt accounts for the rising proflle of the high 
energy fraction, it does not correctly model the behaviour below 20 GeV , instead giving 
a result that is signiflcantly too high. This is not particularly surprising given the rough 
and ready nature of our toy model for the magnetic fleld. 



4.2 An oscillating current sheet 



In the previous section we considered a simplifled magnetic fleld with a flat current sheet. 
This is clearly not realistic - we know that the current sheet oscillates and extends to 



rise to be slower and extend to lower energies. We will present a careful analysis of the proton to anti- 
proton flux in a future paper but for now we note that this model requires that there be no anti-proton 
excess in contrast to many dark matter explanations of PAMELA. 



large heliolatitudes. Even at solar minimuni the current sheet extends up to around 15° 
(from the Wilcox Solar Observatory [17]). The tilt varies considerably from month to 
month. However we will consider an idealised situation where the tilt remains constant 
throughout the heliosphere to qualitatively understand the effect of the tilted current 
sheet on the high energy cosmic ray flux. 



Particle Paths with an Oscillating Current Sheet 
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Figure 4: A schematic diagram of the particle paths in the magnetic field with an oscillating 
current sheet. Particles of both signs can penetrate the center of the solar system along the ecliptic 
plane once their momentum is large enough that their Larmor radius exceeds the wavelength of 
the current sheet. Particles with a large Larmor radius that leave the angular volume described 
by the ecliptic have paths that are equivalent to the scenario of a flat current sheet just as before, 
providing an excess of one charge over the other. 

We show a schematic diagram of high energy particle paths in Fig. |4j In modulation 
theory [TH] the effect of the oscillating current sheet is incorporated taking it to be a fiat 
sheet that gives rise to an overall drift velocity in the plane. The oscillations of the current 
sheet have the effect of increasing the distance travelled by particles traveling along the 
current sheet and consequently reducing the overall drift velocity. At low energies this 
results in a drift velocity that is lower than the solar wind speed and the inward drift of 
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the cosmic rays is countered by the convection of the solar wind. 

This works well as long as Tg is smaller than the wavelength of the current sheet 
(6.5 AU) at all points in the heliosphere. This is true for sub-GeV particles considered in 
classical heliospheric modulation. However as soon as Vg is greater than the wavelength a 
particle can pass directly along the ecliptic - as shown in Fig. |4} This is charge symmetric 
close to the ecliptic. This allows particles of both species to propagate deep into the 
heliosphere. As the Larmor radius goes as 1/r, and the wavelength of the HCS remains 
constant with r, a particle with a given energy will penetrate up to the point where 
Tg becomes less than the wavelength. At this point a positively charged particle (in 
the current solar phase) will continue to travel inwards along the current sheet, and a 
negatively charged particle will not. At a large enough energy, particles of both species 
will travel to the inner heliosphere in this manner. 

The second effect to note is that for high energy particles that travel outside the 
volume containing the current sheet we still have the same asymmetric paths that cause 
positrons move to the center of the solar system and cause electrons to move towards the 
termination shock. This is to be expected, as when r^ is larger than the amplitude of 
oscillation there are many paths on which the positrons cross the current sheet only once 
per rotation. For these particles the oscillation of the current sheet is irrelevant and the 
effect should be identical to the case of a flat current sheet. For these particles we should 
see the effect that we considered earlier where positrons have a direct path to center of 
the solar system and electrons do not. 

To accurately account for the variation in the cosmic ray species due to these complex 
paths will require a detailed computer model of the particle propagation. This is underway 
and will be presented in a future work. For now we construct an approximation to 
highlight the main effect of the oscillating current sheet: 

• High energy particles that arrive at the ecliptic with relatively large incident angles 
will cross the current sheet only once per half rotation - at all points along their 
path of travel. In this case the current sheet is identical to a fiat current sheet. 

• High energy particles that have a gyroradius r^ > 6.5 AU can pass radially through 
the volume containing the oscillating current sheet. This is true for particles of 
both signs and is also true in both directions. For these particles there is no charge 
dependent excess of one species over the other. 

We model this effect by imposing a fiat reduction in the fraction of the particles that 
can be considered to be crossing a fiat current sheet. 

/' = sin0-5, /'>0 (11) 

where we expect 5 to be proportional to sin6'i as the reduction is directly due to the 
latitudinal extent of the current sheet. 
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This does not account for the increase in the local flux that we should see above a few 
GeV in both electrons and positrons from their passage directly through the oscillating 
heliosheet along the ecliptic. However these increases are entirely charge symmetric so 
they will only have a small effect the ratio of fluxes. The more significant effect is to reduce 
the fraction of positrons that can be considered to be crossing a fiat current sheet, and to 
effectively impose a low energy threshold below which no positrons can be considered to 
be traveling along a fiat current sheet. 



Taking this into account we now write the ratio of fluxes to be: 

^iAi/(e+) {A + f)FsoAu{e-^) 



^lAl/(e^ 



{{A + f')FsoAu{e+) + AFsoAuie- 



(12) 



Fit to PAMELA data with an oscillating current sheet 
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Figure 5: We fit the PAMELA data with a model that incorporates a varying current sheet. At 
low energies we do not expect a good fit as conventional heliospheric modulation effects dominate 
here. 



We show the result in Fig. |5} This fit requires the local density to be 15% of the 
interstellar density and for 6 = 0.09. The oscillating current sheet has the effect of 
moving the rise in the positron fraction to larger energies, and decreasing the overall 
magnification effect of the lense. 
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Even though we have used a number of broad approximations the predictions are 
robust. It is true that there is an asymmetry in the propagation paths of particles of 
different signs at these energies. It is also clear that this asymmetry can only become 
significant once the Larmor radius of the particles is large with respect to the latitudinal 
extent of the current sheet, and once it does become significant it rises with energy. We 
only need to assume that there remains a radial gradient in the density of cosmic rays at 
these energieq^to predict a rising positron fraction between 10 and 100 GeV. 

This allows us to make two predictions. Firstly, that the excess will disappear once we 
leave our current solar minimum. Secondly, we note that the tilt of the current sheet has 
a substantial effect on our model as it determines the energy at which the effect switches 
on. The tilt of the current sheet changes over time and as such provides us with an 
opportunity. We would predict a variation of the excess of positrons that correlates with 
large scale variations in the tilt of the current sheet, delayed by the time it takes for those 
variations to propagate out into the wider heliosphere. 



5 Conclusions 



We have proposed that the PAMELA positron excess is not down to annihilating dark 
matter or new astrophysical sources, but instead results from the configuration of the 
heliospheric magnetic field. The ordered nature of the magnetic field on large scales 
creates a lense that allows particles of one sign to free stream into the center of the solar 
system whilst particles off the opposite sign travel out. This effect rises with energy. It 
naturally occurs in the correct energy range for the correct sign of particle to account 
for the observed positron excess. Though our model of the solar system magnetic field 
is rough, we can already make two clear predictions. Firstly, the PAMELA result will 
disappear when we go through solar maximum (with a suitable delay to account for 
the time it takes the conditions at the sun to propagate out through the solar system). 
Secondly, we predict that there will be fluctuations in the observed excess that correlate 
(again with a delay) to the tilt angle of the heliospheric current sheet. 



'*The radial gradient in cosmic ray densities should disappear for energies of a few TeV when particles 
free stream through the heliosphere. It is tempting to use this to explain the Fermi excess in cosmic ray 
electrons. In this case the bump at 500 GeV would be the observation of the unmodulated spectrum 
whereas the low energy data would represent the locally suppressed flux. To fit the Fermi data in this way 
would indicate the local flux at 10 GeV to be 60% of the unmodulated interstellar flux - 4 times the flux 
we use to fit the PAMELA signal. We note that there are many sources of uncertainty in our model that 
could account for this discrepancy and we will require that our more complete description should account 
for this difference. There also remains the possibility that there is an extra primary component to the 
electron and positron spectrums but that such an excess would be smaller than previously considered. 
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